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Abstract 

We propose a simple prescription to calculate the solar neutrino survival 
probability Pee in the quasivacuum oscillation (QVO) regime. Such prescrip- 
tion is obtained by matching perturbative and exact analytical results, which 
effectively take into account the density distribution in the Sun as provided 
by the standard solar model. The resulting analytical recipe for the calcu- 
lation of Pee is shown to reach its highest accuracy (|APee| < 2.6 x 10~^ in 
the whole QVO range) when the familiar prescription of choosing the solar 
density scale parameter rg at the Mikheyev-Smirnov-Wolfenstein (MSW) res- 
onance point is replaced by a new one, namely, when tq is chosen at the point 
of "maximal violation of adiabaticity" (MVA) along the neutrino trajectory 
in the Sun. The MVA prescription admits a smooth transition from the QVO 
regime to the MSW transition one. We discuss in detail the phase acquired 
by neutrinos in the Sun, and show that it might be of relevance for the studies 
of relatively short timescale variations of the fluxes of the solar u lines in the 
future real-time solar neutrino experiments. Finally, we elucidate the role of 
matter effects in the convective zone of the Sun. 
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I. INTRODUCTION 



The solar neutrino problem [Q], emerging as a deficit of the observed solar neutrino rates 
with respect to standard solar model (SSM) predictions [p],|9|,p!0| , can be explained 
through neutrino flavor oscillations |]ll|, possibly affected by the presence of matter [p!^ -[T^ 
(see [0,0 for reviews of oscillation solutions and [T^JT8[] for recent analyses). The analysis 
of neutrino oscillations requires a detailed study of the evolution of the flavor states (z/g, i^x), 
being any linear combination of z/^ and u^-, along the neutrino trajectory. In the simplest 
case, the active states (fg, i^x) are assumed to be superpositions, through a mixing angle 
uj G [0,7r/2], of two vacuum mass eigenstates (z/i,z/2), characterized by a mass squared gap 
m\ > (see, e.g.. 



5m 



TTln 
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The corresponding neutrino evolution equation in the flavor basis. 




H 




involves then the following vacuum (v) and matter (m) terms in the Hamiltonian H, 

H = + Hm 

k ( — cos 2uj sin 2uj 
2 Y sin2ti; cos 2^7 

with the usual definitions for the neutrino wavenumber in vacuum, 

k = 6m^/2E , 

and for the neutrino potential in matter, 

V{x) = V2GFNe{x) , 



(1) 

(2) 
(3) 

(4) 
(5) 



Ne{x) being the electron density at the point of neutrino trajectory in the Sun, located at 
radial distance x from the Sun center .Q The Ne{x) distribution is usually taken from SSM 
calculations P,p!0|. 

The Hamiltonian H is diagonalized in the matter eigenstate basis (z^im, ^2m) through the 
mixing angle in matter Um defined by 



km sm2ujrn = k sin 20; , 
km cos 2ujm = k cos 2uj — V 



(6) 
(7) 



km being the neutrino wavenumber in matter, 



km = k y (cos 2uj — V/ ky + sin 2u! . 



(8) 



The solar neutrinos reaching the Earth move practically radially in the Sun. The effect of off- 
center trajectories is negligible for our results, see Appendix A. 
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In this work, we focus on approximate solutions of Eq. (|l|) for neutrino flavor transitions at 
relatively small values of 5m? / E [pO|- 24 1 , characterizing the so-called quasivacuum oscillation 
(QVO) ^5H2|Jl3|Tg 



regime^, 
QVO 



5rr?/E < 10"^ eVVMeV . 



(9) 



In order to define our goals more precisely, we recall that, both in the QVO regime and in the 
"Mikheyev-Smirnov-Wolfenstein" (MSW) [|12|,|13| regime, which takes place at 5m? ^ few x 
10^^ eV^/MeV, the solar Ve survival probability at the Earth surface can be written as 



where Pq is the average probability. 



Po + Pi 



2~ Pcj COS 2uj COS 2uj^ , 



and Pi is an oscillating term |^ (see also 



Pi 



Pc(l - Pc) cos2uj^ sin2co'cos($2i - $22) 



(10) 



(11) 



(12) 



The main ingredients of the above equations are: {i) The mixing angle in matter uj^ at 
the z/e production point Xq, with uj^ ~ 7r/2 in the QVO regime; {ii) The level crossing (or 
jump) probability Pc of the transition from at x = xq to vi at the surface of the Sun 
(x = Rq)] and the phase 
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22; 



(13) 



accumulated on the v path from Xq to the Earth surface, x 
$2i; ^ = 1;2, have a simple physical interpretation pll 



L = 1 A.U.,0where the phases 
"^•2* = ^2i{xQ,L) = aTg[A2i{xo, L)], 
being the probability amplitude of the transition V2m{xo) — ^ I'iiL) between the initial 
matter-eigenstate i'2m{xo) and the final vacuum- mass eigenstate Vi{L). The phase ^ can also 
be decomposed (see, e.g.. Appendix A) as a sum of a "solar" phase acquired on the path 
in the Sun (x G [xc-R©]) and a "vacuum" phase acquired in vacuum (x G [i?0,L]), 



= is + k{L - Po) 



(14) 
(15) 



^The lower part of the QVO range {5m? /E ^ 5 x 10"^'^ eV^/MeV) corresponds to the vacuum 
oscillation (VO) regime [^. In this work, we do not explicitly distinguish the VO range, but 
simply treat it as part of the QVO range. Let us note that VO solutions of the solar v problem 
are disfavored by the current data [^,17,18|. 



^In this work we consider only solar matter effects. Earth matter effects are not relevant in 
the QVO regime, while in the MSW regime in which effectively Pi ~ they can always be 
implemented through a modification of the expression for P^e — Po, see the discussion in [p7|]. 
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As a consequence of Eqs. ([IOD-(p!5D, solving the neutrino evolution equation (|l]) basi- 
cally reduces to calculating and ^s- This task can be accomplished through numerical 
integration of Eq. (|^) for x G [0:0,-^0], leading to "exact" solutions.^ However, suitable 
approximations to the exact results are also useful, both to speed up the calculations and 
to clarify the inherent physics. 

The starting point of such approximations is usually the analytical solution of Eq. (|I]) 
for the case of exponential density Ne{x) in the Sun p(]|pT[] (see also [PT|-|5^). Deviations of 
the SSM density from the exponential profile can then be incorporated by an appropriate 
prescription for the choice of the scale height parameter 



ro{x) 



d '^'^ 
— \nNe{x) 
ax 



(16) 



characterizing the realistic change of the electron density along the neutrino trajectory in 
the Sun. The exponential density analytical results are derived assuming that tq = const. A 
well-known and widely used prescription for precision calculations of the jump probability 
Pc in the MSW regime is to use the "running" scale height parameter ro{x) in Eq. (|T6|), 
where for given 5m? / E and for uj < 7r/4, x is chosen to coincide with the MSW resonance 
point, X = Xres PSj. 

In the present work, we show how the analytical solution (Sec. II) and the "running 
resonance" prescription for tq (Sec. Ill) can be smoothly extended from the MSW range 
down to the QVO range, and we give a justification for our procedure. The extension is first 
achieved by matching the resonance prescription to a perturbative expression for P^, in the 
limit of small k (Sec. III). A more satisfactory match is then reached (Sec. IV) by replacing 
the resonance prescription, tq = To(a;res), which can be implemented only for u < 7r/4, with 
the "maximal violation of adiabaticity" (MVA) prescription, tq = 'ro(xmva), where for given 
Sw? / E and any given oj G [0,7r/2], Xmva is the point of the neutrino trajectory in the Sun 
at which adiabaticity is maximally violated (or, more precisely, the adiabaticity function 
has a minimum). We also show how the phase can be easily and accurately calculated 
in the QVO range through perturbative expressions at O(fc^), and discuss the conditions 
under which might play a phenomeno logical role (Sec. V). In Sec. VI it is shown that the 
perturbative results essentially probe the low-density, convective zone of the Sun. Our final 
prescriptions for the calculation of Pc and ^s, summarized in Sec. VII, allow to calculate Pee 
using Eqs. (p!0|)-([I^) with an accuracy |APee| ^ 2.6 x 10^^ in the whole QVO range. All 
technical details are confined in Appendixes A-E. 

While this work was being completed, our attention was brought to the interesting work 
|36| , where the QVO range is also investigated from a viewpoint that, although being gen- 



erally different from ours, partially overlaps on the topic of adiabaticity violation. We have 
then inserted appropriate comments at the end of Sec. IV. 



^In this work, high precision numerical results (which we will call "exact") for SSM density |P,1C| 
are obtained through the computer codes developed in p?|]. 
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II. ANALYTICAL FORMS FOR AND 

In this section we recall the analytical expressions for Pc and C,s, valid for an exponential 
density profile, 

iV,(x) = iV,(0)exp(-x/ro) , (17) 

where A^e(O) and ro are derived from a best fit to the SSM profile P,p!0[ (see Fig. 1 in [p7[1), 

iVe(O) = 245 mol/cm^ , (18) 

ro = i?0/lO.54 . (19) 

Assuming Ne{x) as in Eq. ([T7|), the neutrino evolution equation (Q) can be solved exactly 
in terms of confiuent hypergeometric functions [^,0,0,^,0. The associated expression 



for Pf. can be simplified by making a zeroth order expansion of the confiuent hypergeo- 
metric functions in the small parameter 

1^1 = ror(i?0) ^ 0.165 , (20) 

and by using the asymptotic series expansion in the inverse powers of the large parameter 

l^ol =ro\/(a;o) >4.5 X 10^ , (21) 

where Xq ^ 0.25 i?© (bulk of the neutrino production zone). Then one obtains the well- 
known "double exponential" form for [pO[] , 

exp(27rro/c cos^ cu) — 1 



exp(27rro/c) — 1 



(22) 



Concerning the phase ^s, to zeroth order in roK(-Ro) one gets a compact formula |21 
(without using the asymptotic series expansion in |zo|~"), 

= —2 argr(l — c) — argr(a — 1) + argr(a — c) 

-rofc ln[roV^(xo)] + k^R^ - x^) , (23) 

where 

a = 1 + ir^k sin^ uj , (24) 
c = 1 + ir^k . (25) 

Both Eqs. (|^,|^) are valid at any a;, including uj > 7r/4.|^ 

Let us note that the expression (|I0[) for the probability Pee with the average probability 
Po and the oscillating term Pi given by Eqs. (|TI1)-(|T5|), (p^) and (|23|) was shown |pTl , p2| to 
assume the correct form in the VO, MSW transition and "large 5m^" (averaged oscillations) 
regimes. 



^We remark that the restriction cos 2u! > made at the beginning of [^ 21], which is equiv- 
alent to take UJ < it/A for the usual choice 6m'^ > 0, was basicahy functional to obtain Pge < 1/2 
in a certain region of the parameter space (as it was implied by the results of the Homestake ex- 



periment). However, such restriction does not play any role in the derivation of the Eqs. (22) and 
(p^) in [pO|j2l|, although this was not emphasized at the time. This has also been recently noticed 



m 
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III. THE MODIFIED RESONANCE PRESCRIPTION 



In this section, we show how the resonance prescription for calculation of Pc (and Pee), 
valid in the MSW regime and in the first octant of a;, can be modified to obtain accurate 
values for Pc also for — and for uj ^ 7r/4. 

The resonance prescription in the MSW range is based on the following approximations: 
{i) oscillations are assumed (and where shown in |Q) to be averaged out, so that effectively 
Pi ~ and becomes irrelevant^; {ii) Pc is taken from Eq. (|22|) , but with a variable scale 
height parameter tq = To(a;) [Eq. (|16])], with x "running" with the resonance (RES) point 

X = Xr 
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To = ro(Xres) , (26) 



where Xj-es is defined by the resonance condition 

COs2uJm{Xres) = , (27) 

when applicable. In the absence of resonance crossing, it was customary to take Pc = in 



the MSW regime (see, e.g., [§^). In particular, in the MSW analysis of |]S5|, Pc = was 
taken in the second octant of u, where cos2co'm < and Eq. (pT]) is never satisfied. 

The indicated resonance prescription for the calculation of Pc is known to work very 
well over at least three decades in 6m'^/E (~ IQ-^-lO"^ eVVMeV) |3|, with a typical 
accuracy of a few xlO~^ in Pee. However, in the lowest MSW decade [dm? /E ~ 10~^-10~^ 
eV^/MeV), such prescription is not very accurate both in the first octant (where it tends 
to underestimate the effective value of tq) and in the second octant p5| (where Pc is small 



but not exactly zero). It was found numerically in that, for small a relatively small 



and constant value of tq (~ P0/I5.4) provided a better approximation to Pq in the first 
octantQ; this observation has been also discussed and extended to the second octant in ||23 



(where tq ~ P0/I8.4 is used for — 0). Here we derive (and improve) such prescriptions 
by means of perturbation theory. 

The key result is worked out in Appendix B, where we find a perturbative solution of 
the neutrino evolution equation (|l|), in the limit of small values fc, by treating the vacuum 
term Hy as a perturbation of the dominant matter term Hm in the Hamiltonian [Eq. (|^)]. 
The 0{k) solution^ can be expressed in terms of the following dimensionless integral, 

= dp exp (i r dp'V{p')RQ] , (28) 



\ Jp 



^The same conclusion is valid for all other oscillating terms (and their phases) in Pee |21,22|. 



^The numerical results for Pq were obtained in by using the density Ne{x) provided by the 
Bahcall-Ulrich 1988 standard solar model 



^The perturbative expansion can be expressed in terms of dimensionless parameters such as kRQ 
or krQ. However, for simplicity, we use the notation 0{k^), k ^ etc., instead of 0{k'^R^), 
kRQ — > 0, etc. 
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where p = x/Rq, and V{x) is given by Eq. (|^). In particular, the effective value of ro (in 
the limit of small k) is given by 2tt~^Rq Im(J^), namely, 

lim ro = -Rq C dp sin ( df^^plGp iVe(p') R^ , (29) 
/s^o vr JO \Jp J 

independently of uo (i.e., both in the first and in the second octant). 

The above perturbative results show that the asymptotic {k 0) effective value for ro 
depends upon a well-defined integral over the density profile Ne{x). Using the SSM profile 
for Ne{x) we find that 

lim ro = i?©/ 18.9 . (30) 

The same value is obtained through exact numerical calculations. 

The appearance of integrals over the whole density profile indicates that, for small k, the 
behavior of Pc becomes nonlocal, as was also recently noticed in We further elaborate 
upon the issue of nonlocality in Sec. VI, where we show that the 0{k) perturbative results 
are actually dominated by matter effects in the convective zone of the Sun {x/Rq ^ 0.7), 
where the function Nf.{x) resembles a power law rather than an exponential. 

In order to match the usual resonance prescription [ro = ro(Xi.es)] with the value ro = 
i?0/18.9 in the regime of small /c, we observe that, for the SSM density distribution [^,101 it 



is ro(x) = -R0/I8.9 at X = 0.904 i?©. Thus, we are naturally led to the following "modified 
resonance prescription," 

^ f '"o(a;res) if a;res < 0.904 i?© , , . 

^° \ Rq/18.9 otherwise , ^ ^ 

where "otherwise" includes cases with > 7r/4, for which Xres is not defined. Such a simple 
recipe provides a description of Pc which is continuous in the mass-mixing parameters, and 
is reasonably accurate both in the QVO range {5m^/E < 10"^ eVVMeV) and in the lowest 
MSW decade {5m'^/E ~ few x lO-^-lO"^ eVVMeV). 

Figure 1 shows isolines of Pc in the bilogarithmic plane charted by the variables^ 6m} jE e 
[10"^°, 10"'^] eV^/MeV and tan^ uj G [10"^, 10] .0 The sohd lines refer to the exact numerical 
calculation of Pc, while the dotted lines are obtained through the analytical formula for Pc 
[Eq. (p2D], supplemented with the modified resonance prescription [Eq. (0)]. Also shown 
are, in the first octant, isolines of resonance radius for XresZ-R© = 0.6, 0.7, 0.8, and 0.904. 
The maximum difference between the exact (numerical) results and those obtained using 
the analytic expression for Pc? Eq- (^21), amounts to |APc| ~ 7.5 x 10~^, and is typically 
much smaller. Since Pc is not a directly measurable quantity, we propagate the results of 



^In all the figures of this work, we extend the brr? jE interval somewhat beyond the QVO range, 
in order to display the smooth transition to the MSW range. 

^*^The variable tan^ lo was introduced in [^ to chart both octants of the solar v mixing angle 
in logarithmic scale. 



7 



Fig. 1 to probability amplitudes observable at the Earth, namely, to the average probability 
Pq [Eq. (|TT|) ] and to the prefactor of the oscillating term Pi [Eq. (p!2|)].n 

Figure 2 shows isolines of Pq for SSM density, derived numerically (solid lines) and by 
using the analytic expression for Pq and the modified resonance prescription (dotted lines). 
The maximum difference is |APo| — 3.4 x 10^^. Figure 3 shows, analogously, isolines of 
Pi/ cos^. The difference amounts to |A(Pi/ cos^)| ^ 12.5 x 10^^. 

From Figs. 1-3, the modified resonance prescription for tq [Eq. (pi])] emerges as a rea- 
sonable and remarkably simple approximation to the exact results for P^, valid in both the 
MSW and the QVO regimes. However, it does not reproduce the exact behavior of Pc with 
the requisite high precision of few % for 5m? / E ~ 0(10^^) eV^/MeV and tan^u; ~ 0(1). 
This difference can be understood and removed, to a large extent, through the improved 
prescription discussed in the next Section. 



IV. THE PRESCRIPTION OF MAXIMUM VIOLATION OF ADIABATICITY 



In this section we generalize and improve Eq. (pi]) , by replacing the point of resonance 
(xres) with the point where adiabaticity is maximally violated (xmva), more precisely, where 
the adiabaticity function has its absolute minimum on the neutrino trajectory. 

Let us recall that the validity of the resonance prescription tq = 'ro(xres) is based on the 
fact that Pc 7^ in a relatively shorth part of the v trajectory, where the propagation is locally 
nonadiabatic. The resonance condition, however, can be fulfilled only in the first octant of 
uj. The most general condition for nonadiabaticy, as introduced already in the early papers 
T3| , |35| , ^ , |42[| on the subject, has instead no particular restriction in uo |^3|. Such alternative 



condition can be expressed, in the basis {vimi i^2m) relevant for the calculation of Pc, in terms 
of the ratio between the diagonal term (±/i;m/2) and the off-diagonal term {±iduJm/dx) in 
the (traceless) Hamiltonian. More specifically, a transition is nonadiabatic if the ratio 

^ kUx) ^ k^sm^2uj / {cos2u-V{x)/ky y^' 
^^""^ ~ \2duJm{x)/dx\ \dV{x)/dx\[ sin2 2cu j ' ^ 

satisfies the inequality 7(0;) ^ 1 at least in one point of the neutrino trajectory in the Sun. If 
7(x) is large along the whole trajectory, 7(3;) ^ 1, the transition is adiabatic. The minimal 
value of 7(x) identifies the point of "maximum violation of adiabaticity," Xmva? 

7(a;mva) = min7(x) . (33) 

We show in Appendix B, that, along the solar u trajectory, the MVA point is uniquely 
defined, for any u in both octants. In particular, such point can be unambiguously charac- 
terized through the condition^ 
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The phase ^ is separately studied in Sec V. 



^^A handy approximation to the MVA condition, which by-passes the calculation of derivatives 
with a modest price in accuracy, is discussed at the end of Appendix B. 
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^ ^c/x^^"" =0 at a; = x^va • (34) 

In the first octant, in general, Xmva differs from x^es, although one can have Xmva — x^es 
in some limiting cases (see Appendix B for a more general discussion). For instance, as 
Eq. ( |32D indicates, the two points practically coincide, Xmva — ^^res, in the case of nonadiabtic 
transitions at small mixing angles (sin^ 2c(j ^1). Indeed, let us consider for illustration the 
simplified exponential case of Tq = const. In this ("exp") case it is easy to find from Eq. (|32D 
that: 

- <^^a = ro In 1 (l + V^l + 8(1 + tan^ 2^)) . (35) 

Obviously, at small mixing angles (tan2 2ti; ^ 1) we have Xmva — a^res- However, this is no 
longer true for the nonadiabatic transitions at large mixing angles in the QVO regime we 
are interested in. In the latter it follows from Eq. (^5[) , Xmva < a;res- 

We discuss in Appendix B the more realistic case of SSM density. As far as the calculation 
of Pc is concerned, it turns out that Xm^a Xj-es in the limits of small u (or of large k). 
Therefore, the MVA condition smoothly extends the more familiar resonance condition in 
both directions of large mixing and of small k, which are relevant to pass from the MSW to 
the QVO regime. The inequality Xmva < x^es, derived for exponential density, persists in the 
QVO range for the realistic case of SSM density, implying that 

ro(a;mva) > ro{Xres) , (36) 

with ro(x) defined as in Eq. (0). As a conseqence, a difference arises in the value of Pc if 
Xres is replaced by Xmva in the prescription for calculating tq. 
Explicitly, our MVA prescription reads 

^ f '^o(a;mva) if a;mva < 0.904 i^g , , . 

\ Rq/18.9 otherwise . ^ ^ 

Figure 4 shows curves of iso-ro obtained through Eq. (|37[) . The value of tq presents weak 
variations (within a factor of two) in the whole mass-mixing plane and, by construction, 
smoothly reaches the plateau tq = -R0/I8.9 for Srn^ < few x 10"^ eV^/MeV. 

Figures 5, 6, and 7 are analogous to Figs 1, 2 and 3, respectively, modulo the replacement 
of the resonance condition [Eq. (pT])] with the MVA condition [Eq. (^TD]. The MVA pre- 
scription clearly improves the calculation of Pc, Pq, and Pi/ cos^, with an accuracy better 
than a few percent in the whole plane plotted: \APc\ ^ 3.7 x 10"^, |APo| ~ 2.5 x 10"^, 
|A(Pi/cosO| ^3.3 X 10-2. 

We conclude that the analytical formula for Pc [Eq. (^2])], valid for an exponential density, 
can be applied with good accuracy to the case of SSM density, provided that the scale height 
parameter tq is chosen according to the MVA prescription, Eq. (0). 

A final remark is in order. We agree with the author of ||3^ about the fact that, in order to 



understand better the behavior of Pc in the QVO regime, the concept of adiabaticity violation 
on the whole neutrino trajectory is to be preferred to the concept of adiabaticity violation 
at the resonance point. However, we do not share his pessimism about the possibility of 
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using the running value ro(xniva) for accurate calculations of Pc. indeed, Figs. 5-7 just 
demonstrate this possibility. Such pessimism seems to originate from the observation that, 



as k decreases, Pc starts to get nonlocal contributions from points rather far from Xmva 
In our formalism, this behavior shows up in the A; ^ limit [Eq. (^Ol)], where, as mentioned in 
the previous Section, the effective value of tq gets contributions from an extended portion of 
the density profile [see also Sec. VI and Appendix C for further discussions]. Our prescription 
(P7|), however, effectively takes this fact into account, by matching the "local" behavior of 
Tq for large k [tq = ro(xniva)] with the "nonlocal" behavior of Tq at small k [tq = Rq/18.9]. 
In conclusion, the MVA prescription, appropriately modified [Eq. (P7|)] to match the k ^ 
limit, allows a description of Pc which is very accurate in the whole QVO range, and which 
smoothly matches the familiar resonance prescription up in the MSW range. 

V. THE PHASE ACQUIRED IN THE SUN 

In this Section, we discuss the last piece for the calculation of Pee, namely, the solar 
phase ^s- As we will see, this phase can significantly affect the quasivacuum oscillations 
of almost "monochromatic" solar neutrinos (such as those associated to the ^Be and pep 
spectra). Indeed, there might be favorable conditions in which the phase (often negligible 
in current practical calculations) could lead to observable effects and should thus be taken 
into account. 

First, let us observe that, in the QVO range, the size of the solar phase C,s is of 0{kRQ), 
as indicated by Eq. (^) for exponential density [^, and also confirmed through numerical 



calculations for SSM density ^7j. Figure 9 shows, in particular, exact results for the ratio 
^s/kRQ, as a function of Sm^/E, using the SSM density. It appears that, neglecting with 
respect to the vacuum phase is almost comparable to neglect P0 as compared with L. 
Remarkably, there are cases in which corrections of 0{Rq/L) are nonnegligible, e.g., in the 
study of time variations of Pee over short time scales [p^ - ^ , |4^ , induced by the Earth orbit 
eccentricity {e = 0.0167). In fact, the fractional monthly variation of L{t) from aphelion to 
perihelion is 2eL/6 = 1.2 Rq. Real-time experiments aiming to detect time variations of the 
z/ flux in monthly bins might thus test terms of 0{Rq/L) ~ 0{^s/^v), as also emphasized 



at the end of Sec. VI in the work [^ 



Secondly, let us recall that the oscillating term cos^ gets averaged to zero when the 
total phase ^ is very large ^2\. The approximation (cos^) ~ holds in the MSW regime, 



but it becomes increasingly inaccurate (and is eventually not applicable) as k decreases 
down to the QVO regime. In order to understand when ^ starts to be observable (at least 
in principle) one can consider an optimistic situation, namely, an ideal measurement of 
Pee with a real-time detector having perfect energy resolution, and monitoring the flux from 
narrowest solar u spectra (the Be and pep neutrino lines). In this case, the most important — 
and basically unavoidable — source of smearing is the energy integration P^j45|j46|j4^j47[ 



over the u lineshape.|3 It has been shown in (see also [^,|^) that such integrations 



effectively suppresses the oscillating term Pi at the Earth through a damping factor D, 



13 



Smearing over the u production zone is irrelevant in the QVO regime [p^, see also Appendix A. 
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calculable in terms of the v lineshape. Figure 8 shows the factor D for the pep line and 
the for two Be lines, characterized by average energies (£") = 1442, 863.1, and 385.5 keV, 
respectively. The Be and pep lineshapes [having 0(1) keV widths] have been taken from 
Q and respectively. Figure 8 proves that ^ is observable, at least in principle, in the 



whole QV range dra^ / E ^ 10^^ eV^/MeV, as far as the narrowest Be line is considered. Of 
course, the observability of ^ becomes more critical (or even impossible) by increasing the 
initial energy spread (e.g., for continuous v spectra) or by performing measurements (like 
in current experiments) with additional and substantial detection smearing in the energy or 
time domain 

With the above caveats in mind, we set out to describe accurately in the whole QVO 
range. This task is accomplished in Appendixes D and E where, by means of the same 
perturbative method applied earlier to Pc, we obtain the 0{k) and 0{k'^) expressions for ^s, 
respectively. The final perturbative result is 

~ 0.130 {k Rq) + 1.67 X 10"^(A; R^f cos 2uj + 0(P) , (38) 

which is in excellent agreement with the exact result for shown in Fig. 9, in the whole 
QVO range. 

The 0{k) coefficient in Eq. ( |38D is just the real part of the integral in Eq. ( pHj) (see 
Appendix D) and is already sufficient for an accurate description of S,s in the QVO range. 
Remarkably, the right magnitude of this coefficient can also be obtained from the analytical 
expression (|23|), which would give — 0.116 kR^ (see Appendix D). We also keep the small 
Oik"^) term in Eq. (pSf), because it neatly shows that starts to become cu-dependent for 
increasing values of 5m? /E, consistently with the exact numerical results of Fig. 9. 

Figure 10 shows the error one makes on the phase ^s, by using increasingly better ap- 
proximations, for two representative values of tan^tu. The error is given as the absolute 
difference between approximate and exact results for ^g-, in units of a period (27r). The 
lowest possible approximation is simply to neglect Rq with respect to L, so that ^ = kL and 

= kRQ, namely, neutrino oscillations are effectively started at the Sun center ("empty 
Sun"). A better approximation is to start oscillations at the sun edge, ^ = k{L — Rq) and 



= 0, by assuming that the high sun density damps oscillations up to the surface ||12|J49 
("superdense Sun"). Such two approximations (and especially the first) are widely used in 
phenomenological analyses. Figure 10 shows, however, that they can produce a considerable 
phase shift (even larger than a full period for the case of "empty" Sun) in the upper QVO 
range. The fact that the SSM density is neither A^^e = nor A'g = oo is correctly taken into 
account through the 0{k) term in Eq. (pSj), and even more accurately through the full 0{k'^) 
expression for in Eq. (|38|), as evident from Fig. 10. 

The errors estimated in Fig. 10 show that different approximations for can affect 
calculations for future experiments. In particular, the use of the familiar "empty Sun" or 
"superdense Sun" approximations can possibly generate fake phase shifts in time variation 



analyses for real-time detectors sensitive to neutrino lines, such as Borexino pOl or Kam- 



LAND According to the estimates in ||52|, such two experiments might be sensitive to 



monthly-binned seasonal variations for 5m? ^ 5 x 10~^. Figure 10 shows that, in the upper 
part of such sensitivity range, the empty Sun approximation gives (and thus ^) totally out 
of phase, as compared with exact results. Analogously, the superdense Sun approximation 
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induces a phase shift that, although smaller than in the previous case, can still be as large 
as 7r/2 (A^s < 0.2 x 27r ~ 7r/2) in the quoted sensitivity range, and can thus produce a big 
difference [Acos^ ~ 0(1)] in the calculation of Pi. 

It is an unfortunate circumstance, however, that the dominant term in C,s is proportional 
to k. This fact implies that, neglecting ~ O.13Oi?0/L in the total phase ^ = C,v + ^s, is 
basically equivalent to introduce a small bias of the kind 6171^ x (1 + 0.130 i?0/L) = 

dm? X (1 + 6 X 10~^). On one hand, such bias is sufficient to produce a substantial difference 
in cos.^ when 8rr? jE approaches 10~^ eV^/MeV, and is thus observable in principle. On the 
other hand, brr? is not known a priori, but must be derived from the experiments themselves, 
and it will be hardly known with a precision of 0(few x 10~^) for some time. Therefore, 
although may produce a big effect at jixtd values of 5m^, it might be practically obscured 
by uncertainties in the fitted value of . 

In conclusion, we have found a simple and accurate expression for the solar phase ^s, 
to be used in the calculation of the total phase ^ = + k{L — Rq). The solar phase 
C,s produces effects of 0{Rq/L), which can be nonnegligible in high-statistics, real-time 
experiments sensitive to short-time variations associated to neutrino lines, such as Borexino 
and KamLAND, as was also emphasized in [0,^ and more recently in In such 

context, we recommend the use of Eq. ( PB[ ) for precise calculations of Pge at fixed value of 
(5m^, although the observability of certainly represents a formidable challenge. 



VI. PROBING THE CONVECTIVE ZONE OF THE SUN 

In this section we elaborate upon the 0{k) perturbative results discussed previously for 
Pc and (and detailed in Appendix C and D). In particular, we show that they are related 
to the density in the convective zone of the Sun, corresponding to p = x/Rq ^ 0.7. 

The 0{k) results crucially depend upon the quantity J^, which is defined in Eq. ( pHD 
as an integral over an oscillating function, having the density Ne{x) as inner argument. 
Numerical evaluation of such integral for SSM density gives = (13.03 + i 8.32) x 10^^ (see 
Appendix C). It turns out that the largest contribution to comes from the outer regions 
of the Sun, where the integrand oscillates slowly, while in the inner regions the integrand 
oscillates rapidly, with vanishing net contribution to the real and imaginary parts of Irj. 
Numerical inspection shows that the value of Irj is dominated by the p ^ 0.7 range, which 
happens to correspond to the convective zone of the Sun |l|. Therefore, it is sufficient to 
consider such zone to estimate 

In the convective zone, the density profile Ne{x) is better described by a power law rather 
than by an exponential function (see, e.g., [§§]). A good approximation to the SSM density 
for p ^ 0.7 is: 

N^ip) ^ N,il - pf-' , (39) 

with p ~ 2.8 and Np ~ 1.4 mol/cm^. By adopting such expression for the density, and using 
the position z = 1 — p, we get the following expression for J^, 

In= I dz exp{iqz^) , (40) 
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where q = VpR^/p ^ 135 and Vp = \/2GpNp. Once again, we note that the above integrand 
gives a very small contribution outside the convective zone [z > 0.3), since the exponent 
becomes large. Therefore, the upper limit can be shifted from 1 to oo without appreciable 
numerical changes, and with the advantage that the result can be cast in a compact analytical 
form, reproducing the exact SSM numerical result with good accuracy: 



oo 



In — dz exp{iqz^) (41) 

= g-i/Pr(l + p-^)e'^ (42) 
~ (13 + i8.2) X 10"^ . (43) 

Using the above equation and the results of Appendixes C and D, the small-A; limit for 
ro and can be explicited as 

^ ~ - g-i/P r(l + p-') sin(7r/2p) , (44) 

Rq it 

-|- ~ q~'/P r(l + p-') cos(7r/2j9) . (45) 

Such results show that the effective values of ro and of at small k are connected to the 
parameters {p, q) describing the power-law dependence of Ne{x) in the convective zone of the 
Sun [Eq. (pQl)]. Therefore, while for relatively large k neutrino oscillations in matter probe 
the inner "exponential bulk" of the solar density profile, for small k they mainly probe the 
outer, "power-law" zone of convection. 

Finally, let us notice that, in the range 0.7 ^ p ^ 0.9, where the "running" MVA 
prescription ro = ro(xmva) is applicable, the power law in Eq. (|39| ) leads to ro(x)/-R0 = 
(1 — p)/{p — 1), so that ro{x)/RQ ~ 1/18.9 at the matching point p ~ 0.904, consistently 
with the prescription in Eq. 



VII. SUMMARY OF RECIPES FOR CALCULATION OF Pee IN THE QVO 

REGIME 

We summarize our best recipe for calculation of Pee as follows. In the QVO regime, for 
any given value of the mixing angle a; and of the neutrino wavenumber k = 5m? /2E, the Uf, 
survival probability reads 

Pee = Pc COS^ LU + {1 — Pc) Sll? UJ + siu 2uJ\J Pe{l — Pc) COS ^ , 

where 

^ = ^, + k{L- Rq) . 
The value of Pc can be evaluated with high precision by using the analytical form 

^ exp(27rroA; cos^ cij) — 1 
exp(27rroA;) — 1 
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provided that the density scale parameter ro is calculated as 




) if x„,va < 0.904 i?0 , 
i?0/18.9 otherwise , 



where Xmva is the point where adiabaticity is maximally violated along the neutrino trajectory 
in the Sun (see Appendix B). Such expression for smoothly joins the more familiar 
resonance prescription when passing from the QVO to the MSW regime. 

The value of the Sun phase can be calculated with high precision by using the 0{k^) 
perturbative result, valid for 5m? / E ^ 10~^ eV^/MeV, 

is ^ 0.130 (A; i?©) + 1.67 x {k Rr^f cos 2uj . 

We have shown that such phase can play a role in precise calculations of Pee in the QVO 
range. It is not necessary to extend the calculation of in the MSW range, where effectively 
(cos ^ ~ 0) , and is not observable even in principle. Notice that the neutrino production 
point xq does not appear in the calculation of both and in the QVO range. 

Finally, we have checked that the above recipe allows the calculation of Pee with an 
accuracy |APee| ^ 2.6 x 10~^ (and often much better than a percent) in both octants of 
uo for the whole QVO range {5m? ^ 10~^ eV^/MeV). Figure 11 shows, as an example, a 
graphical comparison with the exact results for P^e at the exit from the Sun (x = P0).|3 
It appears at glance that our recipe represents an accurate substitute to exact numerical 
calculations of Pee, in the whole QVO range 5m? / E < 10~® eV^/MeV. 



VIII. CONCLUSIONS 

We have worked out a simple and accurate prescription to calculate the solar neutrino 
survival probability P^e in the quasivacuum oscillation regime. Such prescription adapts the 
known analytical solution for the exponential case to the true density well as to 

the perturbative solution of neutrino evolution equations in the limit of small 5n? jE. The 
accuracy of the prescription is significantly improved (up to 2.6 x 10~^ in Pee in the whole 
QVO range) by replacing the familiar prescription of choosing the scale height parameter tq 
at the MSW resonance point by a new one: rg is chosen at the point of maximal violation of 
adiabaticity (MVA) along the neutrino trajectory in the Sun. Such generalization preserves 
a smooth transition of our prescription from the QVO to the MSW oscillation regime, where 
the two prescriptions practically coincide. We show that at sufficiently small k = 5m'^/2E 
in the QVO regime, the effective value of tq is determined by an integral over the electron 
density distribution in the Sun, Ne{x): 

lim ro = -Rq dp sin ( dp V2Gf NJp') Rq] , 
fc^o vr JO \Jp J 



^^The function Pee at the Earth (x = V) can not be usefully plotted in the range of Fig. 11, due 
to its rapidly oscillating behavior. 
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where p = x/Rq. The main contribution in the above integral is shown to come from the 
region corresponding to the convective zone of the Sun, x ^ 0.7 Rq. Thus, if quasivacuum 
oscillations take place, solar neutrino experiments might provide information about the 
density distribution in the convective zone of the Sun. We also discuss in detail the phase 
acquired by neutrinos in the Sun, whose observability, although possible in principle, poses 
a formidable challenge for future experiments aiming to observe short timescale variations 
of fluxes from solar u lines. 
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APPENDIX A: GENERAL FORM OF Pee AND ITS QV LIMIT 

In this Appendix we recall the derivation of the basic equations for Pee given in the 
Introduction, together with further details relevant for Appendixes B-E. 

Given an initial solar state ^'e = (1)0)^ (T being the transpose), its final survival 
amplitude Aee at the detector can be factorized as 

A^xo, L) = f/Jc, ^e (Al) 

_ f cosu; sin Co" \ / e*^"/^ 
I — sin uj cos / I 




(A2) 



where the U matrices act as follows (from right to left): (i) t/Jo rotates (at x = Xq) the 
initial state Nf, into the matter mass basis; (ii) Um is a generic SU{2) parametrization for 
the evolution of the matter mass eigenstates {I'lm, ^2m) from the origin (x = Xq) up to the 
exit from the Sun {x = Rq) where = Vi, in terms of the so-called crossing probability 
Pc = P{i^2m ^i) and of two generic phases a and (3] (in) evolves the mass eigenstates 
in vacuum along one astronomical distance L, with S,v defined in Eq. ([T5|); and {iv) finally 
rotates the mass basis back to the flavor basis at the detection point {x = L).^ 



^During nighttime, one should insert a fifth matrix to take into account the evolution within 
the Earth, which we do not consider in this work (focussed on solar matter effects in the QVO 
range). Earth matter effects can always be added afterwards as a calculable modification to Pq, 
since they turn out to be nonnegligible only in the MSW regime, when oscillations are averaged 
out and Pi ~ 0, see [^] and references therein. 
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After some algebra, the resulting v^, survival probability can be expressed as the sum of 



an "average term" Pq plus four "oscillating terms" |21,30 



where 



Pi 
P3 



n=0 



+ - Pcj COS 2uj^ COS 2uj , 

COS 2uj^ sin 2uj\J Pdl — Pc) cos(^^ + it — a — (3) 

sin2co'^ cos2uj^JPc{l — Pc) cos(a; — (3) , 

i sin 2ul^ sin 2a;Pc[cos(^ - 2a) + cos(^„ - 2/3)] , 



P4 = - sin 2ci;|^ sin 2uj cos(^^ — 2a). 



(A3) 

(A4) 

(A5) 
(A6) 
(A7) 

(A8) 



It has been shown in pT| , p2| that the last three terms P2,3,4 can, in general, be safely neglected 
for practical purposes, so that one can take 



Pee — Pq + Pi 
which, together with the identification 

= TT - /3 - a 



(A9) 



(AlO) 



leads to Eqs. ([T0|)-(|T2|).Q 

In the QVO regime relevant for our work {k < 10~^ eV^/MeV), the negligibility of P2,3,4 
is evident from the fact that k/V^Xo) ^ 4 x 10~^ for xq < O.25P0, so that the P2,3,4 prefactor 
sin 2uo^ is negligibly small, while the nonvanishing terms Pq^i can be written as: 



,QVO 



pQVO ^ 2 sincj cosco'a/Pc(1 — Pc) cos^ . 



cos^ ujPc + sin^ ujil — Pc 



(All) 
(A12) 



It was shown in ||2^ that the above equations are not spoiled by Earth matter effects,[^ as 
they turn out to be negligible in the QVO range — a fortunate circumstance that considerably 
simplifies the calculations. 



^^One can make contact with the ^>ij phase notation of |21,22,3^] through the following identifica- 
tions: $11 = -a + ^„/2; $12 = (3- ^t,/2; $21 = vr - /? + ^„/2; and $22 = a- ^„/2. The notation in 
the present work expUcitly factorizes out the contribution of the vacuum phase for x £ [Rq , L] . 



^^Equations ( |A11 ) and ( A12| ) in this work coincide with Eq. (28) in |27|, modulo the identification 
Pq = Pc, valid in the QVO regime. 



16 



We have verified the apphcabihty of the approximations in Eq. ( |A11|) and ([A12|) in the 
QVO regime, by checking that our exact results do not change in any appreciable way by 
setting coscj^ = from the start. In particular, the results of the numerical integration of 
the evolution equations vary very little by forcing the initial state to be z/2m rather than Vf,. 
We have also verified that none of our figures changes in a graphically perceptible way in 
the QV range, by moving the point xq within the production zone (xq ^ O.25i?0) not only 
radially but also for off-center trajectories. Therefore, in the following appendixes, we will 
neglect corrections of 0{k/V{xQ)) in the QVO regime, and just set 

xo = , (A13) 
ul ^ 7r/2 , (AM) 

from the start, without any appreciable loss of accuracy. 

Given Eq. (|A^), and the fact that the initial state can be taken, to a high degree of 
accuracy, u^. = z/2m, the probability amplitudes to find the two mass eigenstates at x = Rq 
iiv = 0) simply read 

= /^e^("-^) , (A15) 



^22 = Vl-^ce*°- (A16) 
In terms of flavor state transition amplitudes A^e and A^x at x = Rq, the jump probability 
reads 

pQvo = |^^^|2^ (A17) 

= lAee]"^ COS^ UJ + \Aex\'^ Sin^ CO — 2Re{AeeAl^) Sinu! COSU! , (A18) 

and the phase readsj^ 

^ a,g(^^^ ^*^) (A19) 

= arg[sinu;cosu;(|Aee|^ - l^exH + Re(AeeA*^) cos2u; + i Im(AeeA*^)] . (A20) 

We conclude this Appendix by quoting useful expressions for some dimensionless quan- 
tities which appear in the calculation of Pc and ^: 

L/Rq = 215 , (A21) 



kRQ = 1.76 X 10^ ( I > (A22) 



Sm^/E \ 



eVVMeV / 
^mol/cm'^ 



VRq = 2.69 X 10^ ( — , (A23) 



V/k = 1.53 X 10^' ' — -j—r . (A24) 



Notice that the a^o-independence of implies that no smearing over the production zone 

is necessary, as also emphasized in [p7|]. 

i^Equations ^KWj ) and { \K2^ can also be obtained from Eqs. (13a)~(13d), (15), and (19) in |3C|. 
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APPENDIX B: THE CONDITION OF MAXIMAL VIOLATION OF 

ADIABATICITY 



In this Appendix we characterize the condition of maximum violation of adiabaticity 
along the neutrino trajectory. We also show that the MVA condition reduces to the familiar 
resonance condition in appropriate limits. 

In the matter eigenstate basis, the neutrino evolution equation 





(Bl) 



is governed by the Hamiltonian 



jji / 1 2 i LO„i/k„i \ ('B2') 

~ T V 2 i Urrjkm I j ' ^ ^ 

where tOm = dojm/ dx. 

The evolution is adiabatic when the diagonal term dominates over the off-diagonal term, 
so that the ratio of such elements [the function 7(0;) defined in Eq. (|32|) ] is large. The 
"maximal violation of adiabaticity" along the neutrino trajectory is reached at the point 
X = Xmva where 7(0;) is instead minimized, 

d'y(x) ,^ , 

-T^ = ^ x = x^^,. (B3 
ax 

The above MVA condition can be rewritten as 

d"^ cos 2uJm 



dx"^ 



^ X = Xmva ■ (B4) 



Since cos2ci;m increases monotonically from its value at the production point (~ —1) to its 
vacuum value (cos 2uj) , the above condition characterizes the (flex) point of fastest increase 
for the function cos2u!m{x), which is thus uniquely defined at any u G [0, vr/2]. 
The MVA condition ( |B4| ) can be compared with the usual resonance condition, 

COs2co'm = ^ X = Xres • (B5) 

which can be fulfilled only for u G [0,7r/4]. To understand under which circumstances the 
two conditions practically coincide (xmva — a^res), it is useful to rewrite the MVA condition 
in the following, equivalent form: 

3sin2c<JmCos2c<Jm(V)^ + /csin2co'V = ^ x = Xmva ■ (B6) 



The above condition reduces to Eq. (P5|) in the case linear density {V = 0) (which, however, 
is not applicable to the Sun). It also reduces to Eq. (p5|) for very small mixing (sin 20; — >■ 0), 
which explains why the simple resonance condition has been very frequently used in the 
literature, to replace the slightly more complicated MVA condition. Notice also that, for 
increasing k, Pc becomes nonnegligible only at small mixing angles. Therefore, as far as the 
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calculation of Pc is concerned, the MVA and resonance conditions also merge at large k. In 
all other cases, it turns out that Xrcs > x^^^. for u G [0,7r/4]. For u G [7r/4, 7r/2], Xres is 
simply not defined. 

We conclude by providing a handy approximation to the value of Xmva, which may be 
used to by-pass the calculation of derivatives of V{x), at the price of a small loss of accuracy 
in the calculation of Pee- As discussed before, the MVA condition characterizes the flex 
point of the curve cosa;m(a:), which increases from —1 to cos2u; in the QVO range. We have 
verified that, in the parameter region where Pc is nonnegligible, the fiex point Xmva is close 
to the half-rise point Xi/2 where cos2a;m assumes the average value between its two extrema, 
defined as 

cos2LJm{xi/2) = -(cos 2^7 — 1) . (B7) 
The above definition leads to the condition 

r (X1/2) = k cos 200 + , (B8) 

V 1 + sin uj 

which clearly reduces to the resonance condition in the small uu limit. Using ro(xi/2) as 
a substitute for ro(xniva) in the MVA prescription [Eq. (pT]) ] provides a final accuracy of 
^ 4.5 X 10~^ in Pee in the whole QV regime (only slightly worse than the accuracy quoted 
at the end of Sec. Vll). 



APPENDIX C: EXPANSION OF AT FIRST ORDER IN k 

In this Appendix we study the limit — for Pe, in both cases of exponential and SSM 
density. Such limit gives an asymptotic value for the effective scale parameter rg, relevant 
for the lower part of the QV regime. The positions in Eqs. (|A13|) and ( |A14|) are adopted. 
We also introduce a normalized radial coordinate, 

p = x/Rq . (CI) 

In the case of exponential density, Eq. ( P^D gives, at first order in k, 

Pc = cos^ uj — —r^k siv? 2uj , (C2) 

where tq = Rq/10M [Eq. (^]. 

In the case of SSM density, we demonstrate that Eq. (|C2|) is formally preserved, modulo 
the replacement 

ro ^ r; = P0/I8.9 . (C3) 

The crucial observation is that, for k ^ 10~^ eV^/MeV and SSM density, we have k -C V{x) 
for most of the u trajectory inside the Sun (p < 0.9) before the rapid density drop at the 
border. This suggest a perturbative solution of Eqs. (||) and (0) with the vacuum term 
playing the role of a small perturbation, as compared with the dominant matter term H^. 
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IS 



The perturbative expansion is expected to be accurate when the second term in Eq. 
small, namely, for Zero -C 1. 

At 0th order in H^, the evolution operator from p = to p = 1 is trivial in the flavor 
basis being diagonal). 



ro(l, 0) = exp (^i dp Hmip)^ 



exp 



-1^(1,0) 







exp 
where 

V{p2,pi) = / dpV{p)RQ . 
J pi 

At 1st order in (i.e., in k), one gets an improved evolution operator Ti, 

Ti(l,0) = To(l,0)[l-zA(l,0)] , 
where the correction matrix A reads 

A(l, 0) = T,-\l, 0) dpToil, p) To(p, 0) 
Jo 



'0 



— COS 2uj 

»i 



sin2ti; / dpe 
Jo 



^ \^sin2o^ dpe-'^^^P'^^ cos2c^ 
By applying Ti to the initial flavor state V'e = (1, 0)^, one gets the flavor state at x = Rq, 



(C4) 
(C5) 

(C6) 

(C7) 

(C8) 
(C9) 



Re 



— cos 2lj 

nl 



o) 2^^® I sin2u; j^' dpe'"^^''"^ 



(CIO) 



where the 0{k) correction is explicited. The basic integral, 

I-q = f dpexp[ir]{l, p)] = C + iS , 

•J 

appearing in Eq. ( |C10|) , can be numerically evaluated using the SSM density, 

C = dp cos dp'V{p')RQ = 0.1303 , 
Jo Jp 

S = dp sin dp'V{p')RQ = 0.0832 . 

JO Jp 

Finally, using Eqs. ( |A18|j(JlOD and keeping only 0{k) terms, one gets 

Pc ~ cos^ u — sin 2u! Re^i'ei'*) 



cos uj kRp)S sin 2uj 

2 ° 



(Cll) 

(C12) 
(C13) 

(C14) 
(C15) 
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which reproduces Eq. ( |C2|) , up to the anticipated replacement 



ro->r[, = 2^i?0/7r = /?0/18.9 . (C16) 

We have independently checked the value r'^ = i?0/18.9 by numerical studies of the /c — > 
limit for the quantity {Pc — cos^ u!)/k, with derived exactly rather than perturbatively.0 
We observe that, for the SSM density, the scale parameter vq^x) happens to be equal to 
the value in Eq. Q at x = 0.904 i?©, 

ro(O.9O4i?0) = Rq/18.9 = r'^ . (C17) 

Such relation proves useful to match our analytical and perturbative results around the 
critical zone of the "knee" of the SSM density, which occurs just at a; ^ O.9i?0. 

As a final exercise, we prove that Tq ~ tq in the exponential limit V{x) = V{0)e~^^'^° , as 
it should be. The proof uses the assumptions of large initial density and small final density, 
as discussed for Eqs. ( pUD and (pi]). Then one changes the integration variable from p to 
= V{0)roe~'^^®^^° in the calculation of S [Eq. (|U13| )], and exploits the fact that is very 
small (large) for p = 1 (p = 0): 

^ = - dp sin dp' V{0) Rq e-^'^°/"o (C18) 
Rq vr JO Jp 

dp sin (V{0) ro 6""^®/"°) (C19) 



2 fi 

TT Jo 

ro 2 f°° , , sin tq 



RqTC Jo Rq 



(C20) 



APPENDIX D: EXPANSION OF AT FIRST ORDER IN k 

In this Appendix we study the small-fc limit of the phase ^s, by using the same strategy 
adopted in the previous section for Pc. In particular, let us consider the 0{k) expansion of 
the analytical phase in Eq. (P^D, 

6~fci?0C", (Dl) 

whereQ 

C' = l-^[7i5 + In roV(O)] = 0.116 . (D2) 

Kq 



^'^The author of [p^j3^ found numerically a similar value for Tq (namely, i?0/18.4), by fitting his 
exact results for small k. In this Appendix we have given an explicit derivation of rg, based on 
perturbation theory, and independent on numerical fits. 

'^^'lE — 0.577 is the Euler constant. 
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We show that Eq. (pi|) , vahd for exponential density, remains formally unchanged also for 
the SSM density, modulo the replacement C — >• C, where C = 0.130 is given by Eq. (|C12| ). 
In fact, using Eqs. ( [A20| ) and ( |C10| ), and keeping only the 0{k) terms, we get 

^2Im(AeeA:j/sin2cu (D3) 
= kRQC . (D4) 

Such result is independently confirmed by studying the A; — > limit of the quantity ^s/kRo,, 
with C,s derived from exact (rather than perturbative) calculations. 

Finally, we show that C ~ C" in the limit of exponential density, as it should be. With 
manipulations similar to those used to derive Eq. ( |C20| ) we get: 

C = C dp cos C dp ^(0) ^0 6"'''^®/^° (D5) 
Jo J p 

~ dp cos {y{<d) ro e"''-^®/"") (D6) 

= 1 - dp [l - cos (f(0) ro e^^^®/"")] (D7) 



Rq Jo 

-l-^[7i? + lnroV(0)] = C' (D9) 

where, in the divergent integral of Eq. ( |D8| ), we have kept only the two leading terms [je 
and InroV^(O)] which do not decrease when the upper limit roV{0) becomes large. 



APPENDIX E: EXPANSION OF AT SECOND ORDER IN k 

In this Appendix, we present an 0{k'^) perturbative expression for ^s, which effectively 
accounts for the cu-dependence of the exact results (as evident in the upper part of the QVO 
range in Fig. 9). 

We basically iterate the calculation of the evolution operator T(l, 0) (see Appendix C) 
at second order in the perturbation term H^. We omit the (somewhat lengthy) derivation 
and present the final O(fc^) result, 

= (kRe) Relrj + {kR^f cos 2uj [Im(/^ - + I'^)/2 - ReJ^ ■ Im/^] (El) 
= 0.130 {kRo) + 1.67 X l{]-^{kRQf cos2cj , (E2) 

where = (13.03 + i8.32) x 10"^ jg defined in Eq. ( |CTT| ), while 1'^ and are defined as 

I'r, = dpp exp[iri{l,p)] (E3) 

= (12.45 + ^7.07) X 10-2 , (E4) 

/; = £ dp dp' exp[^ 7^(1, p')] (E5) 

= (0.58 + U. 25) X 10^2 ^ (E6) 
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and have been numerically evaluated for SSM density. The 0{k) term reproduces, of course, 
the results given in the previous section. 

The 0{k'^) behavior of Eq. ( |E2| ) is not captured by the analytical phase written in the 
form of Eq. (^), which contains only odd powers in a fc-expansion. The reason can be 
traced to the fact that Eq. (^) was derived in [^,0 under the assumption of zero density 
at X = -R0 (which is not exactly fulfilled for exponential density). Preliminary studies 
indicate that removal of such assumption gives the correct O(fc^) behavior of the analytical 
phase [including the cos2u; dependence displayed in Eq. (^)]. Such studies are proceeding 
and will be presented after completion. 
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FIG. 1. The crossing probability Pc, calculated for the SSM density in the parameter space 
[brr^jE^ tan^o;). The three lowest decades in brr^jE characterize the quasivacuum oscillation 
regime. Solid curves: exact numerical results. Dotted curves: approximate results, using the 
analytical resonance prescription. Also shown are isolincs of resonance radii. The radius x/Rq = 
0.904 represents the matching point of the resonance prescription with the perturbative results 
valid for small values of Srr? jE. 
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FIG. 2. The average survival probability Po, as obtained through exact numerical calculations 
(solid curves) and through the analytical resonance prescription (dotted curves) . 
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FIG. 3. The oscillating term prefactor Pi/ cos^, as obtained through exact numerical calcula- 
tions (solid curves) and through the analytical resonance prescription (dotted curves). 
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FIG. 4. Isolines of the effective density scale parameter ro (in units of Rq), as derived through 
the prescription of maximum violation of adiabaticity (MVA), matched with the perturbative result 
valid at small values of 5m? /E (tq/Rq — > 1/18.9). 
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FIG. 5. The crossing probability Pc, as obtained through exact numerical calculations (solid 

curves) and through the analytical MVA prescription (dotted curves). Also shown are isolines of 
MVA radii. The radius x/Rq = 0.904 represents the matching point of the MVA prescription with 
the perturbative results valid at small values of 6m? /E. 
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FIG. 6. The average survival probability Po, as obtained through exact numerical calculations 
(solid curves) and through the analytical MYA prescription (dotted curves). 
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FIG. 7. The oscillating term prefactor Pi/ cos^, as obtained through exact numerical calcula- 
tions (solid curves) and through the analytical MVA prescription (dotted curves). 
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Damping foctors for Be ond pep lines 




6mV<E> (eV'/MeV) 



FIG. 8. Damping factors for the oscillating term cos^, as calculated for the Be and pep so- 
lar neutrino lines {{E) being their average energy). The damping factors completely suppress 
oscillations at the Earth for 5m^/ {E) above the quasivacuum range. 
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FIG. 9. The solar phase ^g, in units of fci?©, as obtained from exact numerical calculations for 
some representative values of u. Notice the a;-independent limit for A; — 0. 
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FIG. 10. The absolute error induced by ctpproxiiriEite Ccilculcitions of the solcir phcise , cis 
compared with the exact numerical calculations, in units of 27r. See the text for details. 
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FIG. 11. Application of our final analytical recipe for Pc and to the calculation of PeeiR®) at 
the exit from the Sun (dotted lines), as compared with the corresponding exact numerical results 
(solid lines). 
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